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In this series of eight papers we present the appHcations 
of methods from wavelet analysis to polynomial approxi- 
mations for a number of accelerator physics problems. In 
this part, according to variational approach we obtain a rep- 
resentation for orbital particle motion in storage rings as a 
multiresolution (multiscales) expansion in the base of well- 
localized in phase space wavelet basis. By means of this 
"wavelet microscope" technique we can take into account 
contribution from each scale of resolution. 



1 INTRODUCTION 

This is the first part of our eight presentations in which 
we consider applications of methods from wavelet anal- 
ysis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], which is based on 
our approach to investigation of nonlinear problems - gen- 
eral, with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum, which are 
considered in the framework of local (nonlinear) Fourier 
analysis, or wavelet analysis. Wavelet analysis is a rela- 
tively novel set of mathematical methods, which gives us a 
possibility to work with well-localized bases in functional 
spaces and with the general type of operators (differential, 
integral, pseudodifferential) in such bases. In the parts 1- 
8 we consider applications of wavelet technique to non- 
linear dynamical problems with polynomial type of non- 
linearities. In this part we consider this very useful ap- 
proximation in the case of orbital motion in storage rings. 
Approximation up to octupole terms is only a particular 
case of our general construction for n-poles. Our solutions 
are parametrized by solutions of a number of reduced al- 
gebraical problems one from which is nonlinear with the 
same degree of nonlinearity and the rest are the linear prob- 
lems which correspond to particular method of calculation 
of scalar products of functions from wavelet bases and their 
derivatives. 



2 ORBITAL MOTION IN STORAGE 
RINGS 

We consider as the main example the particle motion in 
storage rings in standard approach, which is based on con- 
sideration in [9]. Starting from Hamiltonian, which de- 
scribed classical dynamics in storage rings H{f,P,t) = 
c{7r^+mQC^}^/^+e(/) and using Serret-Frenet parametriza- 
tion, we have after standard manipulations with truncation 
of power series expansion of square root the following ap- 
proximated (up to octupoles) Hamiltonian for orbital mo- 
tion in machine coordinates: 
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Then we use series expansion of function f{pa) from [9]: 

f{p.) = /(O) + /'(OK + /"(0)p^/2 + ... = p.- 
p'^/{2jq) + . . . and the corresponding expansion of RHS 
of equations corresponding to (1). In the following we take 
into account only an arbitrary polynomial (in terms of dy- 
namical variables) expressions and neglecting all nonpoly- 
nomial types of expressions, i.e. we consider such approxi- 
mations of RHS, which are not more than polynomial func- 
tions in dynamical variables and arbitrary functions of in- 
dependent variable s ("time" in our case, if we consider our 
system of equations as dynamical problem). 

3 POLYNOMIAL DYNAMICS 

The first main part of our consideration is some variational 
approach to this problem, which reduces initial problem to 
the problem of solution of functional equations at the first 
stage and some algebraical problems at the second stage. 
We have the solution in a compactly supported wavelet ba- 
sis. Multiresolution expansion is the second main part of 
our construction. The solution is parameterized by solu- 
tions of two reduced algebraical problems, one is nonlinear 



and the second is some linear problem, which is obtained 
from the method of Connection Coefficients (CC). 

3. 1 Variational Method 

Our problems may be formulated as the systems of ordi- 
nary differential equations 

dx.Jdt ^ fi{xj,t), (i, j = 1, (2) 

with fixed initial conditions Xi{0), where fi are not more 
than polynomial functions of dynamical variables Xj and 
have arbitrary dependence of time. Because of time dila- 
tion we can consider only next time interval: < t < 1. 
Let us consider a set of functions $j (t) = Xidyi/dt + fiUi 
and a set of functionals 

F^ix) = / 'P,it)dt - x,y, \l (3) 
Jo 

where yi{t){yi{0) = 0) are dual variables. It is obvious 
that the initial system and the system 

F^ix) = (4) 

are equivalent. In the following parts we consider an ap- 
proach, which is based on taking into account underlying 
symplectic structure and on more useful and flexible ana- 
lytical approach, related to bilinear structure of initial func- 
tional. Now we consider formal expansions for Xi^yc 

x,{t) = xm + J2 >HVk{t) y^it) = J2vjMt), (5) 

k r 

where because of initial conditions we need only ipk (0) = 
0. Then we have the following reduced algebraical system 
of equations on the set of unknown coefficients A*^ of ex- 
pansions (||): 

^Mfe.A,^-7r(A,) = (6) 

k 

Its coefficients are 

Mfcr = / 'P'kit)fr{t)dt, (7) 

Jo 

ll = / Hx,,t)^r{t)dt. 
JO 

Now, when we solve system and determine unknown 
coefficients from formal expansion we therefore obtain 
the solution of our initial problem. It should be noted if we 
consider only truncated expansion (|]) with N terms then we 
have from ^ the system of iV x 71 algebraical equations and 
the degree of this algebraical system coincides with degree 
of initial differential system. So, we have the solution of 
the initial nonlinear (polynomial) problem in the form 

N 

= 2^,(0) + ^AfXfe(i), (8) 
fe=i 



where coefficients \\ are roots of the corresponding re- 
duced algebraical problem (||). Consequently, we have a 
parametrization of solution of initial problem by solution 
of reduced algebraical problem (Q). The first main problem 
is a problem of computations of coefficients of reduced al- 
gebraical system. As we will see, these problems may be 
explicitly solved in wavelet approach. The obtained so- 
lutions are given in the form (|^), where X}^{t) are basis 
functions and are roots of reduced system of equations. 
In our case Xk [t) are obtained via multiresolution expan- 
sions and represented by compactly supported wavelets and 
A^ are the roots of corresponding general polynomial sys- 
tem (H) with coefficients, which are given by CC construc- 
tion. According to the variational method to give the reduc- 
tion from differential to algebraical system of equations we 
need compute the objects 7^ and fjLji, which are constructed 
from objects: 

<J^ = [ X^{T)dT, V^j = ( X, (r)Xj (r)dT, 

Jo Jo 

^^J^ = J X,'(T)X,(T)dr, (9) 
Pklj = / Xk{T)Xl(T)Xj{T)dT 

Jo 

for the simplest case of Riccati systems (sextupole approx- 
imation), where degree of nonlinearity equals to two. For 
the general case of arbitrary n we have analogous to it- 
erated integrals with the degree of monomials in integrand 
which is one more bigger than degree of initial system. 

3.2 Wavelet Computations 

Now we give construction for computations of objects (9) 
in the wavelet case. We present some details of wavelet 
machinery in part 2. We use compactly supported wavelet 
basis (Fig. 1, for example): orthonormal basis for functions 
in L2(r). 




Figure 1 : Wavelets at different scales and locations 

Let be / : R — > C and the wavelet expansion is 
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If in formulae (10) cjk = for j > J, then f{x) has an 
alternative expansion in terms of dilated scaling functions 
only f{x) = J2eez cji^ji{x). This is a finite wavelet ex- 
pansion, it can be written solely in terms of translated scal- 
ing functions. Also we have the shortest possible support: 
scaling function DN (where N is even integer) will have 
support [0, — 1] and N/2 vanishing moments. There 
exists A > such that DN has AA^ continuous deriva- 
tives; for small N,X > 0.55. To solve our second asso- 
ciated linear problem we need to evaluate derivatives of 
f{x) in terms of fix). Let be = d"iy9^(a;)/dx". We 
consider computation of the wavelet - Galerkin integrals. 
Let /''(x) be d-derivative of function f{x), then we have 
/'''{x) = Tlit'^ifii^)^ ^ii'^ values f^ix) can be expanded 
in terms of (p{x) 



i{x) = ^ X„iip,n{x), 



(11) 



where Am = / (pi{x)(pm{x)dx are wavelet-Galerkin in- 
tegrals. The coefficients A„i are 2-term connection coeffi- 
cients. In general we need to find {di > 0) 

CXD 

Ktt = J ]K(^)dx (12) 



For Riccati case (sextupole) we need to evaluate two and 
three connection coefficients 
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According to CC method [10] we use the next construc- 
tion. When N in scaling equation is a finite even positive 
integer the function 'p{x) has compact support contained in 
[0, A^ - 1] . For a fixed triple {di ,d2,d3) only some Af^^'^^ 
are nonzero: 2~N<e<N-2, 2 - N < m < 
N-2, \£-Tn\ < N-2. There are Af = 3N^-9N + 7 
such pairs {£, m). Let A'^^''^'^^ be an M-vector, whose com- 
ponents are numbers A^^''^''^. Then we have the first re- 
duced algebraical system : A satisfy the system of equa- 
tions {d — di + d2 + d-^) 

ji^^did2d3 _ 2^^d. p^di_d2di (\^) 
P 

By moment equations we have created a system of M + d+ 
1 equations in M unknowns. It has rank M and we can ob- 
tain unique solution by combination of LU decomposition 
and QR algorithm. The second reduced algebraical sys- 
tem gives us the 2-term connection coefficients. For non- 
quadratic case we have analogously additional linear prob- 
lems for objects (12). Solving these linear problems we 
obtain the coefficients of nonlinear algebraical system (6) 



and after that we obtain the coefficients of wavelet expan- 
sion (8). As a result we obtained the explicit time solution 
of our problem in the base of compactly supported wavelets 
with the best possible localization in the phase space, which 
allows us to control contribution from each scale of under- 
lying multiresolution expansions. 

In the following parts we consider extension of this ap- 
proach to the case of (periodic) boundary conditions, the 
case of presence of arbitrary variable coefficients and more 
flexible biorthogonal wavelet approach. 

We are very grateful to M. Cornacchia (SLAC), W. Her- 
rmannsfeldt (SLAC), Mrs. J. Kono (LBL) and M. Laraneta 
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